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This paper examines the thermodynamic stability of chiral domain walls and vortices - topological 
defects which can exist in chiral superconductors. Using London theory it is demonstrated that at 
sufficiently small applied and chiral fields the existence of domain walls and vortices in the sample 
is not favored and the sample's configuration is a single domain. The particular chirality of the 
single-domain configuration is neither favored nor disfavored by the applied field. Increasing the 
field leads to an entry of a domain-wall loop or a vortex into the sample. The formation of a straight 
domain wall is never preferred in equilibrium. Values of the entry (critical) fields for both types of 
defects, as well as the equilibrium size of the domain-wall loop, are calculated. We also consider a 
mesoscopic chiral sample and calculate its zero-field magnetization, susceptibility, and a change in 
the magnetic moment due to a vortex or a domain wall entry. We show that in the case of a soft 
domain wall whose energetics is dominated by the chiral current (and not by the surface tension) 
its behavior in mesoscopic samples is substantially different from that in the bulk case and can be 
used for a controllable transfer of edge excitations. The applicability of these results to Sr2Ru04 - 
a tentative chiral superconductor - is discussed. 



I. INTRODUCTION 

Chiral superconductors belong to an exotic class of 
physical systems whose many-body ground state carries 
nonzero current and hence breaks time-reversal symme- 
try. Because the current can assume two time-reversal 
connected directions, the ground state of a chiral super- 
conductor is doubly degenerate (chiral) . This degeneracy 
opens a possibility for the existence of extended topolog- 
ical defects - domain walls - which connect regions of 
opposite chirality and exist along with conventional de- 
fects such as vortex lines. 

It was recently suggested that graphene at specific dop- 
ing can support chiral superconductivity.^^ Another ten- 
tative candidate for a chiral superconductor is Sr2Ru04 
below 1.5 K (see Refs. |3i4) for review) which is corrob- 
orated by /iSR^l the Kerr effect ,121 and phase-sensitive 
measurements.'^ The candidacy of Sr2Ru04 is, how- 
ever, undermined by the fact that several attempts to 
detect a surface magnetic field generated by the chiral 
current J^^Hm qj. ^-j-^g magnetic moment associated with 
thenP^have not yielded a positive result. 

One of the possible explanations aimed to cut this Gor- 
dian knot of seemingly contradicting observations is to 
assume the presence, on a mesoscopic scale, of an alter- 
nating chiral domain structure which leads to the sub- 
stantial field cancellation. Previous studies reported in 
the literature have focused either on the calculation of 
the domain wall properties such as surface tension and 
accompanying chiral current (see, e.g., Refs. I15I16P or 
on finding the distribution of magnetic fields assuming a 
particular d omain s tructure without attempts to justify 
the latter.ESHHEZIiIl 

In this work we use London theory to address the ques- 
tion of the thermodynamic stability of several domain 
configurations for a simple sample's geometry (such as 
cylindrical), where demagnetizing effects can be easily 
taken into account. We show in particular that in small 
applied fields the equilibrium domain configuration cor- 



responds to a single-domain sample and that the specific 
chirality of this domain is neither favored nor disfavored 
by the applied field. Upon increasing the applied field 
either a domain wall loop or a vortex enters the sample; 
the configuration in which a domain wall forms a straight 
line which terminates at the edges of the sample is never 
favored in equilibrium. Values of the entry (critical) fields 
as well as of the domain wall loop size are given in terms 
of model parameters such as the magnitude of the chiral 
current and a domain wall surface tension. 



Motivated by recent cantilever magnetometry mea- 
surements in Sr2Ru04 (Ref. fHj) . we also consider a meso- 
scopic chiral sample with a hole for which we calculate 
zero-field magnetic moment, susceptibility, and a mag- 
netic moment change due to a vortex or domain wall 
entry. We demonstrate that in the mesoscopic limit the 
size of a loop formed by a very soft domain wall can be 
controlled with the applied field and speculate that such 
an effect can give rise to a controllable transfer of edge 
excitations such as Majorana modes. 



The paper is organized in the following way. In Section 
\n\ we give a phenomenological description of a chiral do- 
main wall used throughout the paper. In Section |III| we 
derive an expression for the Gibbs energy of a chiral su- 
perconductor with topological defects in the strong type 
II limit. In Section |IV| we apply results of the previous 
section to several domain configurations in macroscopic 
samples. This is the core section of the paper. In Sec- 
tion |V] we focus on a mesoscopic geometry. Section [Vl| is 
devoted to overall conclusions. Appendix [K\ contains the 
solution of the London equation for a two-domain circular 
cylinder of arbitrary dimensions with a hole. Appendix 
|B] contains some useful results involving modified Bessel 
functions. 
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II. DOMAIN WALL DESCRIPTION 

Even in the simplest case of an isotropic chiral super- 
conductor the domain wall structure can be quite com- 
plicated and is in general determined by the interplay be- 
tween the material - Ginzburg-Landau - parameters and 
the domain wall geometry. To describe chiral domain 
walls we use a simplified model in which it is modeled as 
a sheet-like object with surface tension a which carries 
chiral current with linear density 1i and is characterized 
by a winding number £dw 

The presence of the current along a domain wall is ne- 
cessitated by the chiral nature of the state in which the 
internal orbital motion of Cooper pairs, while compen- 
sated in the bulk, produces a nonzero charge current i 
on a boundary with vacuum and 1i on a boundary with 
another domain. The presence of such a current leads 
to a discontinuity of magnetic field (or, rather, magnetic 
induction) across the domain wall. 

A general description of a domain wall requires several 
winding numbers, which reflects the multicomponent na- 
ture of the underlying chiral state. The winding number 
relevant to our model, £dw, has a physical meaning of a 
net flux (in units of flux quantum, cf. below) generated 
by the domain wall's chiral and screening currents in an 
infinite superconducting medium.!^ Defined in this way 
the winding number f^w depends on the geometry of the 
domain wall and, in general, need not be an integer. 

The domain wall surface tension a complements our 
treatment of topological defects by specifying its intrinsic 
energy per unit area. Although existing calculations seem 
to indicate that for a chiral p + ip superconductor cr > 0, 
the author is not aware of a general proof that would 
exclude the opposite. 

We consider a model in which both the sample and the 
domain wall are translationally invariant along the direc- 
tion of the applied field and focus only on two domain 
wall configurations - a straight line and a circle - which, 
due to their high symmetry, admit straightforward ana- 
lytical treatment. The main difference between the two 



where $o = hc/2e (< 0), A is the vector potential, ps 
is the superfluid density and 9 is the phase of the order 
parameter. The volume integration extends over the re- 
gion of space occupied by the superconductor and § d^s 
denotes the integration over its surface. The above re- 
sult, derived under the main assumption of the London 
approximation ~ uniform superfluid densitjS^' - is valid 
for a superconductor of an arbitrary geometry and relies 
only on the use of Maxwell's equations and Gauss's theo- 



configurations is that whereas a circular domain wall cre- 
ates a nonzero net flux {£dw 7^ 0) the net flux created by 
a straight domain wall vanishes (^dw = 0). It should also 
be pointed out that both a and i will in general be dif- 
ferent for the two configurations; we do not indicate such 
a difference explicitly, unless otherwise stated. 

We note in passing that the domain structure in chiral 
superconductors need not be similar to that in ferromag- 
netic materials because the latter do not exhibit field 
screening. 



III. SURFACE REPRESENTATION OF THE 
GIBBS POTENTIAL FOR A CHIRAL 
SUPERCONDUCTOR 

Let us start by considering a superconducting sample 
placed in a uniform external magnetic field. The distri- 
bution of currents and fields in the sample is a function 
of the applied field and, in thermal equilibrium, can be 
found through minimization of the corresponding Gibbs 
potential, defined as'^*' 

G = F, + ^y"dV (B2-2B.H), (1) 

where H and B are magnetic field and induction respec- 
tively; the volume integration extends over the space oc- 
cupied by the superconductor and over any cavities con- 
tained in it. The free energy Fs of the sample, which by 
our definition excludes the field energy given by the B^ 
term in eqn. ([T]), may contain terms describing kinetic 
energies of charge and spin currents,'22l spin-orbit inter- 
action energy,!^ effects of kinematic spin polarization,!^ 
etc. 

In London theory the free energy of a superconductor 
is approximated by the kinetic energy of supercurrents 
described by the superfiuid velocity Vg. The sum of the 
kinetic energy of supercurrents and the magnetic field 
energy can be written in the following form!^ 



(2) 
I 

rem. The convenience of such a representation relates to 
the fact that for relevant geometries the surface integra- 
tion is usually more straightforward to perform than the 
volume one. Moreover, using eqn. ^ one can avoid di- 
rect calculation of the magnetic field contribution which 
is usually quite cumbersome. 

Eqn. ^ can also be used in the presence of topologi- 
cal defects if they are treated in the following way: The 
volume integrals should exclude regions of nonuniform 



.(i.,.^W + ^B^)^-i|lf »s . (B X V«, - ^» . (B X A), 



3 



superfluid density associated with the defects, while the 
surface integrals should be complemented by an integra- 
tion over the surface that encloses the excluded volume. 
The "missing" contribution to the free energy can be ac- 
counted for by introducing defects' surface energy F^, 
which can be computed from a more general description 
(e.g., Ginzburg-Landau theory). Such an approach leads 
to the following representation of the free energjP^ 



:-d ^ 



(3) 



where the integral now excludes regions with nonuniform 
superfluid density associated with the defects. In the ex- 
treme type II limit A/^ ^ 1, this approach should give 
a quantitatively good approximation for the Ginzburg- 
Landau energy of a superconductor with topological de- 
fects, whereas in the marginal case A > ^ one might hope 
to get a qualitatively reasonable description. 



A. Gibbs potential of a two-domain cylinder 

For a general sample's geometry, evaluation of the sur- 
face integrals in representation ^ is complicated by the 
spatial dependence of the magnetic field. One exception 
is a geometry that has a translational symmetry along 
the direction of the applied field (i.e., a cylinder with an 
arbitrary cross section). In this geometry the value of 
the field on the sample's surface is constant, which can 
be seen from the application of Ampere's law to a rect- 
angular contour with a side parallel to the field. 

Motivated by recent cantilever magnetometry mea- 
surements on mesoscopic annular Sr2Ru04 samples,l^we 
consider a circular hollow cylinder with an onionlike do- 
main structure shown in Fig. 1. The hole which is char- 
acterized by an integer winding number is provides, for 
small applied fields, the only place where vortices can 
reside,*^ and the two-domain configuration is the sim- 
plest one in which a reduction of the total magnetic mo- 
ment can be achieved (as observed in Ref. H?]). 

Results of this section can also be used for cylindrical 
samples with an arbitrarily shaped cross section, pro- 
vided the distance between the defects and the boundary 
is much larger than A. In this limit, as shown in Sec- 
tion |IV[ a circular domain wall configuration is favored 
energetically over a configuration in which a straight do- 
main wall runs across the sample and terminates on the 
sample's boundaries. 

For a cylindrical geometry with the axis parallel to 
the applied field Hq we have B = H and expression ([!]) 
simplifies to 



G = R 



1 

8^ 



d?r (H2 -2H-H„), 



(4) 



To make use of eqn. ( 4 ) we notice that the onionlike ge- 



ometry shown in Fig. 



consists of three surfaces: 1 - the 




FIG. 1: Onionlike two-domain configuration of a chiral sample 
used in the calculation of the Gibbs potential ([gJi. Arrows 
indicate the direction of the chiral currents. 



and 3 - the outer surface of the sample. The surfaces 
are characterized by their respective radii Rj. For the 
domain chirality shown in Fig. [l] surfaces 1 and 3 carry 
counterclockwise chiral current i while the domain wall 
2 carries a clockwise current 2i. 

Owing to the presence of the chiral current, the mag- 
netic field across each surface is discontinuous. Let us 
denote the field values on the inner (— ) and outer (+) 
sides of the surface j as Hj±. In this notation H^j^ and 
Hi- are equivalent to the applied field Ha and to the field 
in the hole Hh, respectively, and the field discontinuities 
are given by 



i?.3- — Hn 



H2+ - H2- - 2l Hh - Hi+ - I (5) 



where we introduced the field jump i = Ani/c which is 
analogous to a domain wall magnetization used, for ex- 
ample, in Refs. I16|17l 

Let denote the total flux through the area limited 



by the surface j and $i 



$j. Then, using the 



inner surface of the sample, 2 - the domain wall surface. 



definition of the Gibbs potential Q and the surface rep- 
resentation of the free energy (|3| given by eqn. ([2]), the 
Gibbs potential of this configuration is given by 

Sng = 8nU-\<Po\is{Hh-Ha) (6) 

- |$0|4w(if2- -Ha+i)- I($21 - $32) - i?a*3, 

where g and / indicate that energies are taken per unit 
length in the direction of the applied field. As discussed 
in Section [nj i^w is a measure of the fiux carried by a 
domain wall in an infinitely large sample and /o- is the 
surface energy of the domain wall per unit length. 

In the absence of chiral currents and domain walls, i.e., 
when both i and £dw are set to zero, expression ^ coin- 
cides, up to an additive constantl^^lj with that obtained in 
Ref. [23] for a hollow nonchiral cylinder. Notice that one 
should not expect G to be of a simple form G cx M ■ Ha 
where M is the magnetic moment of the cylinder since, 
in general, in the absence of the applied field M 7^ 0. 

The application of eqn. (|6| requires knowledge of fields 
and fluxes in the system in terms of the applied fleld and 
parameters ig, £dw, h and Rj. Such knowledge can be 
obtained from a solution of the London equation. The 
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general solution of the London equation for an onionlike 
geometry is given in Appendix [A] We now proceed to the 
analysis of the "macroscopic" limit of this solution, where 
our results have simple analytical form. In Section |V] 
we relax the "macroscopic" constraint and consider this 
geometry for a sample with arbitrary dimensions. 

IV. STABILITY OF TOPOLOGICAL DEFECTS 
IN MACROSCOPIC LIMIT 

In this section we calculate the Gibbs potential in the 
macroscopic limit when all relevant distances, such as the 
size of the sample, the distance between defects, and the 
sample's boundary are larger than A. In this limit the 
precise shape of the boundary and the defect's location 
relative to it are irrelevant. We consider the four defect 
configurations shown in Fig. [2ja): the Meissner state (no 
defects), the state with a vortex, and the states with 
a domain wall shaped as a loop or as a straight line. 
The first three configurations can be obtained as limiting 
cases of the onionlike geometry of Fig. [T] and hence are 
described by eqn. ([6]). In the macroscopic, limit Gibbs 
potentials of more complicated defect configurations such 
as combinations of those mentioned above can be written 
in a similar way. 

A. Meissner state 



cylinder and can be shown to hold even when the screen- 
ing is geometrically limited as in the mesoscopic settings 
considered in Section IVl 

While it is natural to expect the invariancc under the 
full time-reversal operation, which in the case of the 
Meissner state involves the reversal of both chiral cur- 
rent i and the applied field Ha, the invariancc under the 
reversal of either i or Ha alone (partial time-reversal op- 
eration) as demonstrated by ^ might be considered a 
surprising feature.'^ An intriguing question is whether 
this feature is just a peculiarity of the cylindrical geom- 
etry which possesses translational invariancc along the 
direction of the applied field or has a broader validity. 
While the author does not have a proof of the latter, a 
plausibility argument can be given that suggests that the 
invariancc under the partial time-reversal operation can 
be expected if the sample has a mirror symmetry in the 
plane perpendicular to the field. 

We also note that, as evident from eqn. chiral cur- 
rents give a positive contribution to the electromagnetic 
energy of the system. This statement should also hold for 
samples with dimensions of the order of or smaller than 
A. For such mesoscopic samples positive chiral contribu- 
tion to the electromagnetic energy may become compara- 
ble with the negative condensation energy whose scale is 
set by the thermodynamic critical field. This mechanism 
may hinder the formation of the chiral superconducting 
state and has to be borne in mind when considering the 
possibility of a chiral pairing in very small samples-l^^ 



The Meissner state of a cylinder corresponds to a 
single-domain configuration with no trapped defects. 
Fig. [2ja). Its Gibbs potential is obtained from eqn. ^ 
by setting fa-, £s, ^dw, ^2, and Ri to zero and is given by 
the expression 



87r.gM = {i-Ha)^ 



M, 



(7) 



where $m = 'i'3 is the net flux through the sample gener- 
ated by both chiral and screening currents. In the macro- 
scopic limit for a sample of circumference P, the Meissner 
flux is given by a plausible expression: 



(8) 



which is a direct consequence of the fact that both the 
applied field and the field created by the chiral current 
are screened over a region of thickness A around the outer 
edge of the sample.^^ Given expression Q, the Gibbs 
potential of a solid macroscopic cylinder in the Meissner 
state takes the following form: 



(9) 



This at first sight counterintuitive result implies that a 
chirality of a single-domain Meissner state specified by 
the sign of i is neither favored nor disfavored by an ex- 
ternal field. Although proven in the macroscopic limit, 
this statement is in fact independent of the size of the 



B. Vortex state 

A vortex state corresponds to a single-domain configu- 
ration with a hole with nonzero phase winding is around 
it. Such a configuration is obtained from an onionlike ge- 
ometry by setting ^dw = 0, i?i = ^ ^ A, and then taking 
the limit R2 — )■ Ri- Using eqn. ^ the Gibbs potential 
of the vortex state relative to that of the Meissner state 
is given by the expression 



8^(<?„ - gu) = -|$o|4(-ff« - 2Ha), 



(10) 



where Hy is the value of the magnetic field on the outer 
side of the surface which defines the normal vortex core; 
in the notation used in eqn. (|6]), Hy corresponds to H2+ 
after taking the limit R2 Ri- In deriving eqn. (10) we 



have neglected the vortex core energy and the flux carried 
by it - a step which is well justified in the extreme type 
II limit used here. 



For large applied fields, the energy difference ( 10 ) is 



negative, which means that the Meissner state is ther- 
modynamically unstable. The critical field for the vortex 
entry is determined by the following equation: 



Hci,v — Hy/2. 



(11) 



The field Hy can be found by either solving the London 
equation in the macroscopic limit or by taking appropri- 
ate limits in the results for the onionlike geometry given 
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FIG. 2: (a) Configurations of the defects considered in Section IV The positive direction of the applied field is out of the 
page toward the reader, (b) Values of winding numbers and chiral currents for vortex, loop and straight domain wall (DW) 
configurations. The size of the vortex core is exaggerated for visual purposes. 



in the Appendix [A] In doing so one obtains 



C. State with a domain wall loop 



47rA 



2 ln(2A/C). 



(12) 



Thus, the value of the first critical field for a vortex en- 
try in a chiral superconductor is independent of both 
the magnitude and the sign on the chiral current i 
and, within logarithmic precision, coincides with Hd for 
nonchiral superconductor^ (see, e.g., Ref. [17]). Alter- 
natively, for a fixed chirality, eqn. ( 12 ) demonstrates the 
absence of the field-reversal splitting of Hd.v It should 
be noted, however, that inclusion of the vortex core en- 



ergy in ( 10 ), which tends to be chirality-dependent,'2Slwill 



result in nonzero field- reversal splitting of Hd^v 

One might wonder what happens if the magnitude 
of the chiral current i is such that the magnetic field 
created in its immediate neighborhood is larger than 
Hci,v While this question cannot be answered within the 
macroscopic approximation used in this section, eqn. ( 12 ) 
suggests that, for applied fields smaller than Hci,v, vor- 
tices generated by the chiral current's magnetic field will 
tend to stay away from the bulk "decorating" edges of 
the sample and domain boundaries (if present). 

Another conclusion that can be drawn from the re- 
sults of this section is that, in thermal equilibrium in the 
absence of the applied field, the total flux trapped by a 
chiral cavity, located at distance much larger than A away 
from the sample's boundary, is zero. This follows from 
setting Ha = in expression ( 10 1 and then minimizing it 
with respect to is- 



We now proceed to the configuration in which a do- 
main wall forms a circular loop, i.e., terminates in the 
sample forming an "island" of opposite chirality. The 
Gibbs potential g^wo of such a configuration is obtained 
from eqn. ([6| by setting £s = and taking the limit 
i?i — > 0. This leads to the expression 



STT^dwo 



2i$2 - 



Ha+l) 



(13) 



where fa is the surface energy defined after eqn. 
is the field on the inner side of the domain's boundary, $2 
is the flux through the area limited by it, and <i>3 is the 
total flux through the sample, which includes the screen- 
ing contribution and the flux created by the domain. In 
the macroscopic limit these quantities can be found by 
taking the appropriate limit in the general solution for 
the field distribution, which is given in Appendix In 
this way we obtain 



H2- 

$3 



47ri?A 



1 

"2 



$o|4w-i27ri?A, 

-^wl-J-ol, 



(14) 



where R is the radius of the domain island [= R2 in the 
notation of eqn. ([6)]. The flux $3 consists of the flux <I>m 
generated by the boundary of the sample [cf. eqn. (|8|] 
and of the flux carried by a domain wall loop, ^dw'I'o- 
Plugging these results into eqn. ( 13 ) yields the following 



expression for the Gibbs potential of the circular domain 
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wall configuration: 



AttRX 

2|$0Kdwi?a, 



(15) 



where gu is the Gibbs potential of the Meissner state 
and a is the surface tension of the domain wall, such 
that fa- = 2TrRa. We first minimize the above expression 
with respect to R which yields the equilibrium size of the 
domain island: 



Ro 



mi 



dw I 



(16) 



AnXJi^ + Attu/X 



Evaluated at Ro the difference ( 15 ) is positive for small 
applied fields; upon increasing the field the difference 



( 15 1 becomes negative at some value of Ha which defines 



the critical field for the creation of a domain wall loop. 
In other words, circular domains such as that shown in 
Fig. will be thermodynamically stable only if the 
applied field Ha exceeds a critical value Hd dv " 
by 



defined 



H, 



l,dv 



47ra-/A. 



(17) 



Notice that unlike Hd^v for a vortex entry, eqn. (12 



the critical field for the domain wall loop entry depends 
on the chiral current i but is independent of the flux 
carried by the defect. The reason for the latter is the i'dw 
dependence of the equilibrium domain size as specified 
by eqn. ( 16 ). 

Let H h denote the chirality arrangement of a domain 

wall loop state shown in Fig. [2j [a). Unlike the Meissner 
state discussed in Section flV A[ the energy of this state, 
eqn. (151, is not invariant under reversal of the applied 
field. Equivalently, for a fixed applied field the energies 
of the H h arrangement and of its time-reversal coun- 
terpart — I — (obtained by changing the direction of the 
chiral currents and the sign of ^dw) are different. In par- 
ticular, for a positive applied field the energy of the — I — 

arrangement is larger than that of -I h. Although the 

energy of the former state can be lowered by adding vor- 
tices, it will still be larger than either the H V arrange- 
ment or a pure vortex state and hence cannot correspond 
to a true equilibrium. 

Putting together eqns. (12 16 and ITtI) leads to the 



following expression for the size of the domain island: 



Ro/X (X Hci,v/Hci,dv 



(18) 



i.e., Ro scales as the ratio of critical fields of the vortex 
and domain entries. Strictly speaking, the formulation 
that leads to this scaling is valid only if i?o/A 1. How- 
ever, because of the exponential falloff of the screening 
currents, one might expect that it is qualitatively correct 
even in the limiting case of relatively small domains when 
Ro/X > 1. 



To estimate the actual value of i?ci,dw for a given ma- 
terial a knowledge of chiral current i and domain wall 
surface tension a is required. These can be calculated 
using Ginzburg-Landau theory and turn out to depend 
on various material and geometrical parameters such as 
Ginzburg-Landau expansion coefficients and the orien- 
tation of the domain wall relative to the crystal axes 
(see, e.g., Refs. I15|16l) . Ignoring for simplicity material 
anisotropy, one can conclude that 



Ana/X = ei 



l$n 



47r2A3^' 



£2 



47r2A4' 



(19) 



where £1^2 are dimensionless parameters w hich in a weak- 
coupling BCS limit are of the order of l.EIESl We now 
define the following parameter: 



>f<j = 4w/(?A) =£iA/(£2e). 



(20) 



As can be seen from eqns. ([16| and ( |17[ ) this parame- 
ter determines whether the energetics of the domain wall 
is dominated by the chiral current ("soft" domain wall, 
Xd 1) or by the surface tension ("hard" domain wall, 
3> 1). In the weak-coupling limit when £1^2 ^ 1, one 
generally expects that Xd w A/^. Provided the weak- 
coupling limit is applicable for Sr2Ru04 (A/^ ^ 1 and 
hence Xd ^ 1), one would expect that Hd^v ^ ^ci.dw 
and the domain size Ro ~ A. However, given the uncon- 
ventional nature of superconductivity in Sr2Ru04, the 
applicability of the weak-coupling results to this mate- 
rial remains an open question. 



D. State with a straight domain wall that 
terminates at the edges 

A straight domain wall configuration that terminates 
at the edges of the sample is qualitatively different from 
the closed configuration discussed earlier. Unlike the lat- 
ter, the total flux carried by a straight domain wall is 
zero (see, e.g., Ref. IT9]) . which substantially changes its 
energetics. Let <I>-|- denote the total flux carried by ± 
domains shown in Fig. 2]ja). The Gibbs potential of such 
a state can be found along the lines that led to eqn. (|6|, 
and is given by 



dw 



I =8^/^-1- Z ($+-$_ )-i/a ($++*-). (21) 



In the macroscopic limit, fluxes "I>± can be easily com- 
puted, which leads to the following expression for g^^^ : 



87r(5dw| - 9m) = 2i?A (? + Ancr/X), 



(22) 



where is the Gibbs potential of the Meissner state and 
R is the length of the domain wall segment. Provided 
that the surface energy cr > 0, the field-independent 
expression (22 1 is always positive, which implies that a 



straight domain wall configuration is thermodynamically 
unstable relative either to the Meissner state or to the 
state with a domain wall loop. 
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FIG. 3: (a) The spatial profile of the Gibbs potential of a mesoscopic chiral cylinder with a domain wall loop as a function 
of its size for several applied fields. The potential is given relative to that of the Meissner state. Two metastable equilibria 
are clearly visible at zero applied field. (b,c) Critical field ffci,dw and the equilibrium size R'^'^ of the domain wall loop at 
Ha = /fcl, functions of the parameter >Cd, defined by eqn, 

the variation of either ei or £2 



(201 



It is assumed that the variation of Xd is entirely due to 
Surprisingly, even for a mesoscopic sample, values of Hci,dv, and R!^^ shown in (b) and (c) agree 



quantitatively well with analytical results (171 and (16 1 obtained for a macroscopic sample. 



E. Summary of Section |IV 



V. MAGNETIC RESPONSE OF A 
MESOSCOPIC CHIRAL SAMPLE 



To summarize, in this section we obtained critical fields 
for a domain wall i?ci,dw, eqn. (17), and a vortex i?ci,-uj 
eqn. (12), entries into a chiral superconductor in the 



macroscopic limit (all relevant dimensions much larger 
than A). It was shown that the preferred domain wall 
configuration is that of a loop whose equilibrium size is 
given by eqn. (16). These results imply, in particular, 
that, at fields above -ffci,dw, a cross section of the do- 
main structure of a macroscopic sample in the direction 
perpendicular to the field is that of a plum pudding - a 
single domain populated by domain islands of opposite 
chirality. 

It was also pointed out that a statement, often encoun- 
tered in the literature, that cooling in the field should re- 
duce domains by biasing the system to one chirality (field 
training) does not refer to the thermodynamic equilib- 
rium. This can be seen from the eqn. (|9| for the Gibbs 
potential of a single-domain sample that does not contain 
terms linear in the chiral current and hence cannot differ- 
entiate between domains of opposite chirality. Upon in- 
creasing the field a topological defect which corresponds 
to the minimal of the two fields Hd.v and i?ci,dw enters 
the sample and for large fields both vortices and domain 
wall loops will be present. 

However, even in thermal equilibrium, the direction 
of the applied field can affect relative chiralities; for ex- 
ample, for a positive field the H h domain wall loop 

configuration of Fig. [2ja) is favored over its time reversal 



In this section we consider the thermodynamic stability 
of domain walls and vortices in mesoscopic chiral samples 
whose relevant dimensions are comparable to X.^^ This 
problem is motivated by recent cantilever magnetometry 
measurements done on small Sr2Ru04 particles.'^^ Al- 
though the main aim of Ref. il4i was to probe the ex- 
istence of half-quantum vortices, it is interesting to ex- 
amine whether the observations reported there shed any 
light on the question of the chiral nature of Sr2Ru04. 

As in Section |IV| we will make use of expression ^ 
to evaluate the Gibbs potential for the onionlike geome- 
try shown in Fig. [T] It is assumed that vortices present 
in the system reside only in the cylinder's hole, which 
guarantees a contour-independent definition of the vor- 
tex winding number ig. This assumption excludes the 
possibility of wall vortices and limits our consideration 
to relatively small applied fields and moderate chiral cur- 
rents (cf. discussion at the end of Section IV B). Notice, 



however, that, as demonstrated in Ref. [141 in a confined 
geometry with geometrically reduced screening the field 
required for the wall vortex entry can be substantially 
larger than the bulk Hci,v [given by eqn. (12)]. 



It is convenient to introducing the notation 

ajk = KoU)Ioik)-Io{j)Ko{k), 
b,k^K2{])h{k)-h{j)K^{k), 



(23) 



where /„ and Kn are modified Bessel functions of rt-th 
order and Ko{j) = Ko{Rj/X), etc.; Ri, R3, and i?2 stand 
for the radii of the inner and outer surfaces and for the 
radius of the circular domain wall loop, respectively (see 
Fig.[l]). 
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The magnetic moment of the cyHnder is given by the 
foUowing expression: 



(24) 



where L is the height of the cyUnder and $3 is the total 
flux through the area limited by outer boundary 3. Un- 
like the macroscopic limit, the values of fields and fluxes 
required to evaluate quantities of interest are no longer 
given by simple analytical expressions and are relegated 
to Appendix [Aj where the general solution of the London 
equation for this geometry is presented. 

In the zero applied field the ground state of the system 
is obtained through minimization of ([6]) and corresponds 
to ^dw = 0, £s = i.e. to a single-domain defect-free state. 
The zero-field magnetic moment Mq of a single-domain 
chiral cylinder is then equal to 

Afo = (iL/At:){ttRI ~ 2^AV&i3) + (25) 

where Xm is the magnetic susceptibility, Xm = OM/dHa- 
The magnetic susceptibility is determined by the sys- 
tem's dimensions i?i and R3 and does not depend on 
the chiral current i: 



AttL ^Xn 



ai3 



4A4 



1 



^1^3 ^13 



(26) 



It is interesting to note that the result ( 26 ) also holds for 
nonzero and £dw, i-e., in the presence of either (hole) 
vortices or a domain wall loop as long as the radius of 
the latter is field independent. In mesoscopic settings the 
independence of equilibrium Ri on Ha can be expected 
for a "hard" domain wall whose energetics is dominated 
by the surface tension a. In the opposite limit of a "soft" 
domain wall whose behavior is dominated by the chiral 
current i and not by the surface tension a one can expect 
significant variations of Ri with if^, as demonstrated 
below. Such variations lead to a deviation of the response 



from the simple linear form described by ( 26 1 



Upon increasing the applied field the system under- 
goes a transition into a state in which either la or ^dw is 
nonzero. The ordering of these events can be determined 
from the comparison of the critical fields required for the 
entry of the defects. The critical field for a hole vortex 
entry is given by 



Hcl,v 



ai3 



l^ol 

27ri?2 bi3-2X^/Rl 



(27) 



Setting i?i — >■ ^ and i?3 — )• 00 in the expression above, 
one recovers the bulk limit given by eqn. (12). 



To find the critical field for a circular domain wall en- 
try one first needs to know its equilibrium size, which 
can be found through the minimization of the Gibbs po- 
tential Q with respect to i?2. Unlike the macroscopic 
limit (Section IV C), the spatial profile of the Gibbs po- 
tential in geometries with constrained screening can be 
quite complicated and may include several metastable 
equilibria [see Fig. [3^], which obstructs transparent ana- 
lytical treatment. Numerical results for -ffci.dw and .Rj'^ 



are given in Fig. ^h) and (c) , where they are plotted as 
a function of the parameter xr^ which characterizes the 
interplay between chiral currents and the surface tension, 
eqn. (20). It has also been checked that the dependencies 
shown inn Fig. [sjjb) and (c) also describe a cylinder with 
a hole, provided that R2^ is constrained to lie between 
Ri and R^ and Hd^dw is constrained by the values of 
which correspond to R2^ = i?i and R"^^ = R3. 

It is important to emphasize that both Hd^v and 
Hci.dw discussed above are computed for a defect-free 
sample. Only one of these fields has a physical meaning; 
for example, if it turns out that Hd.v < -ffci,dw then the 
value of the latter needs to be recalculated in the presence 
of a vortex. 

The entry of a defect into the sample leads to a jump 
in the magnetic moment. Such a jump can be evaluated 
using the results of Appendix |A] and is given by 



AM,, 



An 



AL 1 



for a hole vortex entry and 



AM 



dw 



l^ol^ 
47r 



dw 



2A^ 

i??6l3 



bi2 
bi3 



(28) 



(29) 



for a circular domain wall entry. While AM^ is inde- 
pendent of chiral current, AMdw depends on the domain 
wall size R2, which is determined through the energy 
minimization and hence implicitly depends on i. In the 
limit i?2 -J> Ri we have AM„/AMdw = AiJAld^ and if 
R2 R3 then AMdw ^ 0. 

We now come back to the case of an extremely soft 
domain wall mentioned earlier. Fig. |4] shows the radius 
of a circular domain wall with cr as a function of 
the applied field. As the applied field is increased the 
domain wall moves continuously from the outer to the 
inner surface of the cylinder.'^ Given the possibility that 
chiral boundaries can carry topologically nontrivial exci- 
tations such as Majorana modes (see, e.g., Ref. 1^5]) one 
may speculate that such a process can be used to perform 
a controllable transfer of excitations between the edges 
of the sample. 



A. Application to Sr2Ru04 

We now turn to the question of the interpretation of the 
results of Jang et al^ in terms of possible chiral super- 
conductivity. Jang et al. reported cantilever magnetome- 
try measurements of a mesoscopic Sr2Ru04 particle with 
approximate dimensions Ri — 390 nm, R3 = 850 nni, and 
L = 350 nm; the magnetic moment sensitivity was of the 
order of 10~^^ e.m.u.'*'^ The range of fields used in the 
measurements was such as to cover the first expected en- 
try for a hole vortex, eqn. ( 27 ) . 
The quantities x?i 



AMy and i?ci,u computed earlier 
are independent of chiral current and the only quantity 
that can be used to estimate i independently of a is Mq. 
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FIG. 4: Radius of a soft domain wall loop as a function of 
the applied field. The demonstrated dependence of R2 on Ha 
might be useful for performing controllable transfer of edge 
excitations. 



Within the noise resolution, measurements reported in 
Ref. _14j did not observe zero- field moment, w hich sets 
the limit Mo < 10""'^^ e.m.u. Using expressions (251 and 



(261 and the dimensions of the sample quoted above, we 
obtain the following upper bound for the magnitude of 
the chiral current: 



i < 10 ^ X i„ 



(30) 



where the weak-coupling value iwc for Sr2Ru04 is ob- 
tained from eqns. (19) by setting €2 — I, X — 200 nm 
and is approximately equal to 1.9 x 10^^ statamp/cm. 



Limit ( 30 ) is consistent with scanning SQUID microscopy 
measurement^^^l where, assuming the domain size of the 
order of 1 fj,, measured i was estimated to be less than 
0.1% of the weak-coupling value. 



It was shown that a preferred configuration of the do- 
main wall is that of a loop; a straight domain wall is 
never favored in thermodynamic equilibrium. Domain 
wall loops can exist in the superconducting bulk only if 
the applied field is larger than i?ci,dw, which depends on 
both the magnitude of the chiral current and the surface 
tension of the domain wall. The critical field required for 
a bulk vortex entry is not affected by the presence of the 
chiral currents. 

We have also considered magnetic response and de- 
fect stability in mcsoscopic chiral samples. It was shown 
that for a very soft domain wall its size can be controlled 
by the applied field. This phenomenon can potentially 
provide a mechanism for a controlled transfer of edge ex- 
citations such as Majorana modes. 

There are several possible extensions to this work 
which can be treated in the general framework outlined 
in Section [mi An obvious one is to generalize the results 
presented here for sample geometries that are not trans- 
lationally invariant along the applied field. In particular, 
it is interesting to inquire whether the conclusion that 
the chirality of a single-domain sample is not favored by 
the applied field holds for other types of geometries. One 
geometry which seems to be analytically tractable is that 
of the Pearl limit, in which the thickness of the sample is 
smaller than the penetration depth. This might be par- 
ticular relevant in connections with recen t sp eculations 
about chiral superconductivity in graphene^. 

One can also consider domain structures that are not 
translationally invariant along the direction of the field. 
Such a possibility which was suggested in Ref. [12] is very 
attractive because it has the potential to reconcile the ab- 
sence of the chiral field in Sr2Ru04 as observed by scan- 
ning measurements with nonzero Kerr and /iSR signals.!^ 

Another possible extension would deal with the inter- 
action between domain wall loops and domain wall loops 
and vortices. The interaction between such defects may 
have measurable signatures in the magnetization curves 
of macroscopic chiral samples. 



VI. CONCLUSIONS 



Let us review the quantitative results of this paper. We 
have considered thermodynamic stability of two types of 
topological defects - vortices and domain walls - which 
can exist in a chiral superconductor. Using the London 
theory it was shown that in the zero applied field a macro- 
scopic chiral sample is either defect-free or has defects 
which are expelled toward the edges. The first situation 
is realized if the chiral currents are small and the second 
requires them to be sufficiently large such that the magni- 
tude of the field created in the immediate neighborhood 
of the edge chiral current is larger than a critical field 
required for a defect entry. 
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Appendix A: Solution of the London equation for an 
onion-like geometry 

In this appendix we consider the solution of the London 
equation for the onionlike circular two-domain configura- 
tion shown in Fig. [l] a single-domain configuration can 
be obtained as a limiting case by moving the domain wall 
to the inner or outer boundary of the sample. 

As mentioned earlier, in the extreme type II limit, the 
magnetic properties of a chiral domain wall can be mod- 
eled by replacing it with a sheet current i if the domain 
wall is on the surface and 2i if it is in the bulk. The cur- 
rent carried by e.g. a surface domain wall is iL, where 
L is the length of the domain wall across the direction 
of the current. In a cylindrical geometry with an exter- 
nal field parallel to the cylinder's axis, a sheet current i 
is equivalent to a boundary condition for the magnetic 
field, 



47r . 

H — H^. = — i . 

c 



(Al) 



where the arrow in the subscripts indicates the direction 
of the current and a dot indicates the side at which the 
field is taken. In the London limit the calculation of the 
magnetic response of a sample with a given domain struc- 
ture is thus reduced to solving the London equation in 
each domain and then matching solutions using appro- 
priate boundary conditions. Let R2 denote the radius of 
the domain wall and i?i and i?3 be the inner and outer 
radii of the cylinder. Using eqn. ( Al I, the boundary con- 



ditions for the field can be written in terms of the fields 
on the domain boundaries: 

i/3_ = Ha+l, H2+-H2- = 2l H,+ = Hn-l (A2) 

where iJg- = i^(i?3 - 0), = H{Ri + 0), = 

H{R2 ± 0) and Ha and Hh are the applied field and the 
field in the hole respectively. In cylindrical coordinates, 
the solution of the London equation can be written in 
terms of Bessel functions /q and Xq: 

r e (i?i,i?2) : H{r) - I \) + c'^^K^ir / X) , 

r e (i?2, i?3) : H{r) = cMr/X) + 4;Ko{r/ X), 

(A3) 

where the constants Cjk and c^j. are determined by fields 
on the domain boundaries. 



C12 = a^i{H2-Ko{l) - Hi+Ko{2)), 

c',2 = a^2{Hi+h{2) - 
C23 = a23i(H3_i^o(2) - H2+Ko{3)), 

C23=«27(^^2+/o(3)-i/3-/o(2)), 



(A4) 



with Qjk defined as 



ajk = Ko{j)Io{k)-Io{j)Ko{k), 



(A5) 



where Ko{j) = Ko{Rj), etc. 



The equations above determine the field distribution 
through yet-unknown values of the fields on the bound- 
aries. To find the latter, one can use additional con- 
straints such as those provided by the Feynman-Onsager 
(FO) quantization condition, which is obtained from the 
London form of the Ginzburg-Landau equation for the 
current: 



is = 



87r2A2 



27r 



(A6) 



where is a phase of the superconducting order parame- 
ter and A is a vector potential of the total magnetic field. 
Using the symmetry of the problem and integrating the 
expression above along a circular contour R, one obtains 



_c_ dH_ 
Att dr 



c|"I>n 



87r2i?A2 



(£+l/|$o| ^ A-dZ), (A7) 



where the current density has been expressed in terms of 
the field derivative with the help of the Maxwell equation. 
The constant I characterizes the order parameter phase 
winding around the integration contour. 

In applying the FO relation (A7), one needs to bear 



in mind that the domain wall, being a phase defect, may 
possess a nonzero vorticity and, hence, along with vor- 
tices, contributes to the winding number ^ for appropri- 
ate integration contours. For a circular contour the do- 
main wall vorticity is denoted as ^dw For a nonzero -^dw 
not only the magnetic field but also the screening current 
experience a jump across the domain wall. This follows 
from writing down (A7) for inner and outer boundaries 
of the domain wall. Taking into account that the flux 
is a continuous function of the integration contour, one 
obtains 



{dH/dr)\^,~{dH/dr)\^_ = 



l$n 



27ri?9A2 



(A8) 



where R2 is the radius of the circular domain wall. 

Recalling that due to the presence of the chiral cur- 
rent the magnetic field itself experiences a jump across 
the domain wall and using ( AS I and ( A3 1 we obtain an 



expression that relates field values on boundaries 1,2 and 
3: 



H2-ai3 = Hi+a2z + -ff3-ai2 - « ^20^12(^23 - 023) 



27rA: 



r^dw 1120237 



where R = R/X and 623 is defined bj^^ 

b^k = K2U)Ioik)~l2{j)K„{k). 



(A9) 



(AlO) 



As a useful check of various identities one can consider 
limiting cases of domain wall 2 moving to either the inner 
or the outer surface of the cylinder: 2 — 1 or 2 — > 3. For 
example in the limit 2 — > 1 relation (A9) becomes H2- = 
Hi^ which implies that in a very tl 
is uniform. In the limit 2 



;hin domain the field 
3 we have H2- = H^- — 2i 
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or using eqns. (A2), i^2^ 
conclusion.*'^ 



, thus reaching the same 



Another equation which relates and H2- can be 
found by writing the FO quantization condition ( A7 1 for 



boundary 1. Taking into account that the flux through 
the hole is determined by the field Hh, which is related 
to as Hi^ = H}^ — i, one obtains 



H 



1+ 



ai2 ^-o . , 2 ^ai2 .. . 

I ^^'s + 7^^^^2- - , (All 

012 7ri?f Rjbi2 012 



where is is the number of vortices trapped in the hole. 
Taking the 2 — > 1 limit in the above equation yields 
Hi^ = H2- confirming the expectation that in a van- 
ishingly thin domain the field is uniform. 

Combining eqns. (A9) and (All) allows one to deter- 
mine the field values in the hole Hh and on the domain 
wall, H2- , in terms of the applied field Ha ■ The field in 
the hole is given by the following expression: 



ai3 \<^o\ 


f «23 
J, 

Ol3 




2 


bi3 ttRI 


7:Rl ^ 


Rlbi3 


— R^ai2 


0l3 


- Rlbi2) ^ 


2 


Rlbi3 



Ha 



(1 - ^2^23), 
(A12) 



and the limiting values of the field-independent constant 
Xi are equal to 



3 : Xi 
1 : Xi 



1 



1 



(ai2 



Rl 



(ai 



R\' 



(A13) 



Since in the limits given above the domain structure de- 
generates to a single domain with positive (2 — > 1) or 
negative (2 — 3) chirality the constant Xi changes sign 
as expected. 

The field on the inner side of the domain wall is given 

by 



H2-=- 

^2 = - 



a23 



01 



613 TTRi 

0-23 
bl3 



^12023 |$0 



5i2 



^13 



27rA2 



Ha + 1X2 



{l-Rlb,2)+^^{l~Rib23). 

bi3 



(A14) 



In the limiting cases the field-independent constant X2 
reduces to 



2^3 
2 ^ 1 



X2 = -1 
X2 = -2 



7-(ai3 + ^). 
O13 itf 



(A15) 



The first of these equations implies, in particular, that 
for a zero applied field Ha — and zero winding num- 
ber is = the field at the inner side of the external 
boundary in the limit 2 ^ 3 is equal to —i, as ex- 
pected in this single-domain zero-field limit. The limit 



2 — 1, however, deserves a more careful consideration. 
In this limit H2- should be related to the field in the 
hole as H2- = Hi^ — i and for the zero-applied field zero- 
winding number case the field in the hole can be read 
from eqn. (A13). Comparing the second of eqns. (A13) 



to the second of eqns. ( A15), one indeed recovers such a 
result. 

Having expressed values of the field on all domain 
boundaries in terms of the applied field, the chiral cur- 
rent, and the winding numbers [eqns. (A12| and (A14|] 



we thus find the field distribution in the entire system 
[eqns. (|A3|, (fAl and ([AS])]. 



To write the Gibbs potential of the system requires 
knowledge of fluxes through areas limited by the domain 
boundaries. Although they can be found by the direct 



integration of ( A3 ) it is convenient to find fluxes sur- 
rounded by domain boundaries using the FO quantiza- 
tion condition. Doing so for the second domain boundary 
at r = i?2 + we get 



023 1^0 I 

27rA2 



(4+4w+$2/|$0|) -J^3-- 



rr -^2 



{b2 



-0,23), 

(A16) 

which can be used to find $2- Alternatively, flux $2 can 
also be expressed in terms of fields H2- and Hi+ which 
corresponds to writing FO for the inner part of domain 
wall 2 i.e. for r = R2 — Q. 
Similarly, for $3 we obtain 



a23|^o| 
27rA2 



(4+4w+*3/|*0|) - -H2++H3- ^ (632+023). 

(A17) 

Notice that fiux $3, i.e. the total flux through the area 
limited by the outer boundary of the cylinder, determines 
its magnetic moment M. 

Finally, flux $1 through the opening is determined by 
the field in the hole: 



$1 = nRiHh. 



1. Single-domain limit 



(A18) 



The results for a two-domain cylinder given earlier can 
be used to find the field values in a single-domain limit. 
This is achieved by moving the domain wall located at 
R2 either to the outer R2 — )■ R3 or to the inner R2 — >■ Ri 
boundary of the sample. Doing the former, one recovers a 
single-domain sample with a negative chirality, whereas 
in the latter case a positive chirality sample with 
extra flux quanta trapped in the hole is obtained. 

Taking the limit R2 — >■ R3 the Gibbs potential of a 
single-domain chiral cylinder with a hole is given by the 
expression 

87r.gs.d = -|$o|4(i^ft-i^a)-i^a*3 + «($3-*i), (A19) 

where to agree with the previous notation the total flux 
through the cylinder is denoted as $3 and $1 is the flux 
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through the hole. The field in the hole is given by 



and 



ai3 |$o| 



2A^ 



-H„ 



6i3 T^Ri i?f&l3 

+ ^(-l + (ai3 + 2AVi?2)/^^3), (A20) 



and for the flux through the area limited by the outer 
boundary we obtain 

(A21) 

Notice that it is possible to obtain the single-domain 
results quoted above directly from the results of Ref. |24| 
for a nonchiral cylinder using the following replacements: 
is- TrRfi/\^(il Ha + i, and Hh+ 1 



R 



^{bkiakj — bkjaki) 



(B2) 



Below we give limiting values of several parameters 
required to obtain the results of Section [V] from general 
formulas given in Appendix [A} 

For i?3 ^ 1 the following hold 



ai3 Kq{1) 



8i?3 



0{Rf)) (B3) 



Appendix B: Some useful expressions involving . n\ (^ , ^ ,^o^D-2^^ it3a\ 

modified Bessel functions ^13 ' ^/2^j,^ SKs ^ ' 

In this appendix we collect some useful results which 
involve modified Bessel functions. The unit of length is 

set to A. gfl3 

Given definitions (23 1 of aij and bij one can show that 631 — > —Ko{l) (l — — — — h 0{R^ )) (B5) 



V2^^ 8i?3 



Ko{Ri)h{R2) + Io{Ri)Ki{R2) = ^(621 + 012) (Bl) 
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Our Id-w is n+ in the notation of Ref. 1191 
Several alternative definitions of the Gibbs potential can 
be found in the literature. They differ from the one used 
here by either subtracting the field energy in the absence 
of a superconductor —-^ J H„ and extending the volume 
integration to the whole space or by including the field en- 
ergy of the space occupied by the superconductor ^ J^^ 
into the definition of Fs- 

It is straightforward to generalize this result for situations 
in which the superfluid density is a steplike function of 
coordinates. 

We drop the "potential energy" term V^d-i/;!) from the con- 
sideration since in our linear theory contribution of such 
term away from defects reduces to a field- and current- 
independent constant. 

The possibility of a half-quantum vortex state for which £s 
is a half-integer can be straightforwardly incorporated in 
the formalism but is ignored here for the sake of simplicity. 
Cf. the footnote on different definitions of G. 
Expression ([s]) for the net Meissner flux holds for both solid 
and hollow samples; in the latter case it is required that 
the distance between the hole and the external boundary 
is much larger than A. 

Such invariance also holds for a state with a straight do- 
main wall discussed in Section IIVDI but is violated for a 
state with a domain wall loop, Section [IV C| 
We use convention in which the electron's charge e is neg- 



ative and the phase 9 of the order parameter is defined as 
i/j = \tp\ exp +iO. This implies that a gauge-invariant phase 
is given by the following combination VS-|-27rA/|$o| where 
A is the vector potential such that for a positive applied 
flux the equilibrium winding number is is negative. 
Assuming that the chiral superconductivity is possible 
in mesoscopic samples. Cf. discussion at the end of Sec- 
tion |IVH 

Notice that depending on the actual value of Hd.v a 
higher-field part of the dependence shown in Fig. [4] may 
correspond to a metastable state. 

This is unprecedented moment sensitivity. Compare with 
the sensitivity of torque magnetization measurements in 
cuprates which is about 10~^e.m.u., see e.g. L. Lu et ai, 
Phys. Rev. B 81, 054510 (2010). 

Kerr measurements reported in Ref. [7] are confined to the 
surface since the penetration depth of the EM radiation 
used in Ref. T is of the order of 100 nm. 
Utilizing well-known relations between Bessel functions it 
is possible to use their alternative combinations e.g. those 
involving Bessel functions of zeroth and first order. In par- 
ticular 623 - 123 = ^{lo{3)Ki{2) + Ko{3)Ii{2)) . Our def- 
initions follow those of Ref. 1241 

Notice that in the limit Rj — >■ Rt: a^k — > and hjk — > 
2/Rj. In general ajk ~ —a-kj but bjk cannot be expressed 
in terms of bkj- 



